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The interaction of a Josephson fluxon with an external magnetic field-induced 
potential in a long Josephson junction is investigated experimentally. The 
thermal activation of the fluxon from a potential well is observed and exper- 
iments probing its predicted quantum properties are discussed. A method for 
engineering a magnetic double-well potential for a fluxon is proposed and the 
use of the coupled fluxon states for quantum computation is suggested. 
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1. INTRODUCTION 

Long Josephson junctions are Icnown as tlie most clean systems to ex- 
perimentally study basic properties of Josephson vortices, often also called 
fluxons.EI The interaction of a Josephson vortex with a potential induced 
by an externally applied magnetic field has been discussed in literature in 
various contexts .H'S Recently, first measurements reporting the thermally ac- 
tivated escape of the vortex from a magnetic field-induced potential well at 
4.2 K have been reported. □ Experiments probing the quantum properties of 
a Josephson vortex trapped in a field-engineered potential are currently in 
progress. 

Our interest in this subject is related to a possibility of using macro- 
scopic quantum states of Josephson vortices for quantum computatiorHH. 
We suggest that two distinct states of a fluxon trapped in a field-controlled 
double- well potential inside a narrow long junction can be used for designing 
a qubit, the basic unit of information in a quantum computer. We propose 
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that by varying the external magnetic field and junction shape one is able to 
form an arbitrary-shaped potential for a fiuxon in a junction; the amplitude 
of this potential can be varied in experiment by tuning the magnetic field. 

In this paper we present our experimental results on the thermal escape 
of a fiuxon from a potential well measured at different temperatures. We 
also discuss the quantum tunneling of a fiuxon and some preliminary ideas 
on the engineering of a double- well potential for a fiuxon. 



1.1. A Fiuxon in an External Potential 

We consider a long Josephson junction of arbitrary shape and constant 
width w. w is smaller than the Josephson penetration depth Aj the 
junction can be considered as quasi-one dimensional. Its electrodynamics 
is described by a wave equation for the superconducting phase difference 
4>{x,i) across the junction 



sm0 = a 



c/-,--7 + /(x,t). (1) 



In this representation the time and the space coordinates are normalized 
with respect to the inverse p asma frequency t = ioj~^ and the Josephson 
length X = xXj, respectivelylil In these units, the characteristic phase ve- 
locity of electromagnetic waves in the junction is unity and is usually called 
the Swihart velocity.El The left hand side of Eq. (||) is the sine-Gordon equa- 
tion (SGE). The right hand side of Eq. (|l|) contains the perturbation terms. 
The term 7 = I/Ic is the bias current normalized to the critical current, 
the damping term due to the quasiparticle resistance p is proportional to 
a. The term f{x,t) describes the interaction of the Josephson vortex with 
a spatial inhomogeneity region in the junction. In absence of perturbations 
the Josephson vortex is described by the exact soliton solution of the left 
hand side of Eq. (||) 



(x,t) =4 arctan 



exp 



X — q{t) 



(2) 



where q{t) is the coordinate of the vortex and u is its velocity^ 

Using perturbation theory to solve the perturbed SGEQ (|i]) one can 
write down an equation of motion for the center-of-mass coordinate q{t) of 
the vortex in the non-relativistic limit (u <C 1) 

mfq + aq+^^^ =0. (3) 
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Because of the magnetic flux quantization, the trapped vortex behaves as a 
topologically stable, particle-like object of mass mj moving under the action 
of external forces. The force acting on the vortex can be expressed in terms 
of the potential U. The bias current through the junction acts as a driving 
force (of the Lorentz type) on the vortex. The pinning force may have 
different physical origins. A pinning center for a Josephson vortex in a long 
junction can be realized by introducing a micro-shorti or a micro-resistori 
in the junction barrier. Pinni ng may also occur due to the interaction of the 
vortex with the junction leadal3 or interaction with parasitic magnetic flux 
trapped in the superconducting films. Alternatively, a vortex can be pinned 
by its magnetic dipole interaction with an external magnetic field.a'U In the 
following, we consider the latter pinning induced by an external magnetic 
field, because this method allows to control the pinning potential for the 
vortex during the experiment. 



1.2. Magnetic Field Interaction 

The interaction of a Josephson vortex with an external magnetic field 
applied in the plane of an annular junction can be taken into account by the 
perturbation terrrolHI 

/(x) = /iAsin(^^) (4) 

in Eq. (|l|), where I = lirr/Xj is the normalized circumference of the Joseph- 
son junction, h is the magnetic field normalized to Hq = $o/(2vr/ior(i), where 
d is the magnetic thickness of the junction, and A is a geometry dependent 
magnetic field coupling coefficient .Hil 

Equivalently, in Eq. (^) the fluxon potential U depends on the bias 
current and the magnetic field as given by 

U = - J (-f(f) + A(j)s:n ■ dx, (5) 

where n is the vector normal to the junction and lying in the junction plane, 
and h is the normalized field vector. In the limit of local interaction (Aj — > 
0), the magnetic field part of the potential can be expressed asU = —nofl-H, 
where p is the magnetic moment of the vortex normal to the boundary of 
the junction. 

In the annular Josephson junction shown in Fig. |l|, the fluxon has the 
lowest energy in the location where fl is directed along H and the highest 
energy in the diametrically opposite location where p is opposite to H. The 
pinning potential for the vortex becomes U{q) ~ —hcos{2iTq/£). Subject 
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Fig. 1. Annular Josephson junction with a trapped Josephson vortex 
(fluxon). Under the action of a bias current I > Ic the vortex moves with 
a velocity u generating a voltage V (x u across the junction electrodes. In 
the static state the vortex is pinned in the direction parallel to the external 
magnetic field as indicated by the dashed arrow. 



to a driving force mediated by the bias current /, the potential is tilted 
proportionally to /, effectively lowering the potential barrier for the fluxon. 
Thus, the fluxon dynamics in an annular junction placed in a magnetic field 
is analogous to the motion of a particle in a tilted wash-board potential and 
can be mapped to the dynamics of the phase in a small Josephson junction.113 



2. THERMAL ACTIVATION 

At high enough temperatures, a Josephson vortex trapped in a potential 
well can escape from the well by a thermally activated process. The thermal 
activation rate can be expressed as 

rth = ^exp(^^-^j, (6) 

where loq is the normalized characteristic small amplitude oscillation fre- 
quency of the vortex at the bottom of the potential well. The potential 
barrier height in normalized units Uq{I) and the attempt frequency ooq{I) 
are calculated from Eq. (^). 8£'o = 8$o jc A jw/(27r) is the rest energy of the 
vortex, governing the energy scale of the problem. 

The activation of a vortex from a potential well can be determined ex- 
perimentally bj^^neasuring the switching current distribution of the Joseph- 
son junction .113113 In this type of measurement, the probability of the vortex 
escape from a pinned state (no dc voltage drop across the junction) to a prop- 
agating state (dc voltage drop proportional to the average vortex velocity) 
is determined in dependence on the bias current. The experimentally ob- 
served critical current distribution P{I) is to be compared with a theoretical 



Annular Long Josephson Junctions in a Magnetic Field 




30 40 50 60 
voltage, V [|j,V] 



90 



Fig. 2. Current-voltage characteristic of a fluxon in the annular Josephson 
junction. The modulation of critical current Ic with the external magnetic 
field H is shown in the inset. 



prediction 

r(/')d/'^, (7) 

which is calculated from the bias current sweep rate dl/dt and the depen- 
dence of the activation rate T on the bias current /. If the escape of the 
vortex is due to a thermally activated process, i.e. when ki,T > hujQUjp, the 
critical current distribution P{I) is determined by Eq. (|6[). 

The thermally activated fluxon escape has been experimentally investi- 
gated in a number of samples. Here, we present data on an annular Josephson 
junction of radius r = 50 fim and width w = 3 fim; its geometry is depicted 
in Fig. ||. The junction has a critical current density of jc ~ 100 Acm~^ and, 
consequently, a Josephson length of A j ~ 30 fim. A single Josephson vortex 
could be trapped in the junction by applying a small bias current through 
the junction during the transition from the normal to the superconducting 
state. The single-vortex state is identified by the lowest quantized resonant 
voltage step (saturating at a about ^ = 84 fiV for this particular sample) 
observed on the current-voltage characteristics of the junction, see Fig. |2|. 
Also, a characteristic critical current modulation with magnetic field, accom- 
panied by a suppression of the critical current by a factor of more than 100 
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Fig. 3. Measured histogram (solid squares) of the critical current distribution 
at T = 6bath K and H = 0.25 Oe and calculated probability distribution P{I) 
(solid line). Experimental normalized activation rate r(I) (open squares) 
and fit (dashed line). 



at zero field, as shown in Fig. y and reported earlier in Ref. 17, is observed 
when a single vortex is trapped in the junction. 

A typical measured switching current distribution P{I) at H = 0.25 Oe 
and T = 6K is shown in Fig. |3[ The distribution is normalized such that 
/ P{I)dI = 1. By inverting Eq. (0) the activation rate r(/) can be de- 
termined from the experimental data. Using a cubic approximation of the 
potential Eq. (^), the expression (In (27rr/a;o))^^^ is expected to be linear in 
the bias current / as indeed observed in experiment, see data (open squares) 
in Fig. ^ From a linear fit to (In (27rr/wo))^^^, both the effective temper- 
ature Tesc of the thermal escape and the critical current in the absence of 
fluctuations Iq can be determined self-consistently. The results of this anal- 
ysis are shown in Figs. ^ and b for the bath temperatures T^ath = 2, 4, 6 K 
and magnetic fields in the range of = . . . 0.5 Oe. The error bars in 
Fig. ^ indicate the uncertainty in the value of T^sc calculated using the 
fitting procedure explained above. 

At high temperatures, the experimentally determined values of Tesc are 
within 1 K the thermal bath temperatures T^ath for all values of the mag- 
netic field. For the lowest temperature, the fitted escape temperature is no- 
tably higher than the bath temperature, which is due to the limited current 
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Fig. 4. (a) Escape temperatures Tgsc (open symbols) and (b) fluctuation-free 
critical current Iq (closed symbols) in dependence on the externally applied 
magnetic field at temperatures T = 2 (triangles), 4 (squares), and 6 (circles) 
K. The respective bath temperatures are indicated by dashed horizontal 
lines. 

resolution of approximately 0.5 in the 1 mA range of the measurement 
setup used for data acquisition in this particular experiment. Moreover, a 
small but systematic increase of the escape temperature for small magnetic 
fields is observed. This effect may be explained by considering a localized 
and magnetic field independent contribution to the pinning potential^, which 
leads to an enhanced activation rate at low fields. In Fig. the magnetic 
field dependence of the fluctuation-free critical currents Iq{H) is plotted. 
The currents are observed to be larger than the critical currents measured 
in the Ic{H) patterns, as expected. 

3. QUANTUM TUNNELING 

In the limit of T — > 0, the thermal activation of the vortex from the 
potential well is exponentially suppressed and eventually the temperature 
independent contribution due to quantum tunneling of the vortex through 
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the potential barrier will be the dominating escape process. By measuring 
the switching current distribution we expect to observe a crossover from 
thermal activation to quantum escape at low temperatures. Experiments 
aiming to observe this effect are currently in preparation. 

The theory of macroscopic quantum tunneling for a single Josephson 
vortex has been recently developed and discussed in detail in Refs. 18,19|. 



Essentially, the quantum tunneling of the vortex can be described in the 
Caldeira- Leggett approached, where the escape rate is given by 

T^^ = Aexp{-B) (8) 

with 

/ D \ 1/2 

A = VoOoJoLUpi—j (l + 0(a)), (9) 
B = 1^(1 + 1.74a + 0(a2)). (10) 

This analysis is valid for an expansion of the potential U{q) up to third or- 
der in the vortex coordinate q. Our numerical and analytical calculations of 
the expected switching current distributions (^) using Eqs. show that 
a cross-over from thermal activation to quantum tunneling can be expected 
at temperatures T « 20 mK for junctions which are several micron wide. In 
these calculations we took into account the effects of dissipation which can 
be modeled using the effective damping parameter a in Eq. ([To|). We con- 
sidered both the quasiparticle tunneling across the junction barrier and the 
quasiparticle impedance of the electrodes as possible origins of dissipation 
(also see Ref. |l8|). At 4.2 K the damping parameter a is typically 10~^ and 
decreases exponentially with temperature. Therefore, this damping mecha- 



nism does not substantially decrease the cross-over temperature. In Ref. 19, 
the dissipation induced by the interaction of the vortex with plasma waves is 
studied, and the resulting decrease of the tunneling rate is found to be small. 
So far no other interaction mechanisms of the system with the environment 
have been taken into account. 

Our calculations show that the cross-over temperature can be increased 
close to 100 mK, by reducing the junction width down to sub-micron size and 
optimizing the junction electrical properties. Most experiments on fluxon 
dynamics conducted so far have been done using long Josephson junctions 
of width w > 3fim. This limit is imposed by the standard photolithographic 
preparation procedure. Recently, we succeeded to fabricate high quality 
ultra-narrow long Josephson junctions based on niobium- aluminum-oxide 
trilayer technology of width down to less than 0.3 /im.Ell In this procedure. 
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Fig. 5. (a) Photograph of a 0.3 fim wide heart-shaped junction fabricated in 
Nb-A10a:-Nb technology, (b) Sketch of heart-shaped junction with magnetic 
field applied along the symmetry axis of the junction. Two energetically 
identical vortex states are indicated by arrows in positions and 1. 

we are using electron-beam lithography for junction definition and cross- 
linked PMMA for junction insulation.^ We have successfully prepared and 
characterized narrow long junctions of different geometries. A photograph 
of a 0.3 fim wide heart-shaped long Josephson junction is shown in Fig. ^a. 
First experiments on the activation of fluxons from magnetic field-induced 
potentials in 0.3 fiui wide annular junctions have been recently performed 
and will be presented elsewhere. 



4. POTENTIAL ENGINEERING 

Since the interaction energy of the vortex with the external magnetic 
field is determined by the shape of the junction according to (^), it is possible 
to design an experiment with, in principle, arbitrarily shaped fluxon poten- 
tial by simply varying the junction geometry. Since the major goal in our 
experiments with quantum fluxons will be a demonstration of macroscopic 
quantum coherence, one needs to realize a double-well potential for a vortex 
in the junction. We suggest to use a heart-shaped Josephson junction, as 
shown in Fig. |5|a. The magnetic field H applied along the symmetry axis 
in the junction plane forms two equal potential wells for the fluxon at loca- 
tions "0" and "1". At these two positions the fluxon magnetic moment /2 is 
directed along H, see Fig ^3. The barrier height between the wells can be 
conveniently controlled by the field amplitude and its angle with the sym- 
metry axis of the junction. The corresponding potential profile at h = 1, 
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Fig. 6. Normalized magnetic field interaction potential U{q) at /i = 1 calcu- 
lated according to Eq. (P). 

calculated from Eq. ^ with the normalized junction length L/Xj ~ 35, is 
plotted in Fig. ^ in units of Eq. We note that there are constraints on the 
potentials that can be generated by changing the junction shape, because 
the shape dependent part of the potential is effectively averaged over the 
vortex size (oc Aj), as can be seen from Eq. (^). In the limit of Aj ^ oo all 
features of the potential disappear. 

In the quantum limit, the two distinct fluxon states in the double well 
potential may be employed as a degenerate two-state system. Under suffi- 
cient decoupling from the environment and provided that the temperature 
and dissipation in the junction are low enough, the superposition of the 
macroscopically distinct quantum states "0" and "1" is expected to be ob- 
served. The preparation of the initial state for a quantum coherence mea- 
surement may be achieved by turning the magnetic field in the plane of the 
junction, and the final state of the vortex can be read out by performing an 
escape measurement from one of the potential wells. 



5. CONCLUSIONS 

We have considered a vortex in a long Josephson junction interact- 
ing with an externally applied magnetic field. Experiments indicating the 
thermal activation of the vortex from the field-induced potential have been 
successfully performed. Moreover, the prospects of observing quantum tun- 
neling of a single vortex from such a potential are discussed as well as possi- 
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bilities of potential engineering in order to form a double-well potential for 
use in a qubit. 
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